submitted to Astro. Part. Phys. 
preprint VAND-TH-99-12 



Signatures for a Cosmic Flux 
of Magnetic Monopoles 



Stuart D. Wick^, Thomas W. Kephart^, Thomas J. Weiler^ 
and Peter L. Biermann^ 

^ Department of Physics & Astronomy, Vanderbilt University, 
Nashville, TN 37235 
^ Max Planck Institute for Radioastronomy, 
D- 530 10 Bonn, Germany 



Any early Universe phase transition occurring after inflation has tire potential to populate the 
Universe with relic magnetic monopoles. Observations of galactic magnetic fields, as well as observa- 
tions matched with models for extragalactic magnetic fields, lead to the conclusion that monopoles 
of mass ~ 10^^ GeV are accelerated in these fields to relativistic velocities. We explore the possible 
signatures of a cosmic flux of relativistic monopoles impinging on the earth. The electromagnetically- 
induced signatures of monopoles are reliable. The hadronically-induced signatures are highly model- 
dependent. Among our flndings are (i) the electromagnetic energy losses of monopoles continuously 
initiate a protracted shower of small intensity; (ii) monopoles may traverse the earth's diameter, 
making them a probe of the earth's interior structure; (iii) in addition to the direct monopole 
Cherenkov signal presently employed, a very attractive search strategy for monopoles is detection 
of their coherent radio- Cherenkov signal produced by the charge-excess in the e"^-e~ shower — in 
fact, Cherenkov-detectors have the potential to discover a monopole flux (or limit it) several orders 
of magnitude below the theoretical Parker limit of 10~^^/cm^/s/sr; (iv) it is conceivable (but not 
compelling) that bound states of colored monopoles may be the primary particles initiating the air 
showers observed above the GZK cutofl'. 



I. INTRODUCTION 

Any breaking of a semisimple gauge symmetry which occurs after inflation and which will leave unbroken a [7(1) 
symmetry group may produce an observable abundance of magnetic monopoles. The initial monopole number density 
can be estimated from the type and temperature of the symmetry breaking phase transition [1]. From this, one easily 
derives the present day flux of monopoles [2]. The monopole mass is expected to be ~ a^^ times the scale Ksb 
of the symmetry breaking. It is noteworthy that the inferred strength and coherence size of existing extragalactic 
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magnetic fields suggest that any monopole with a mass near or less than 10^'* GeV would have been accelerated 
in magnetic fields to relativisitic velocities. On striking matter, such as the earth's atmosphere, these relativistic 
monopoles will generate a particle cascade. It is the purpose of this paper to calculate the shower signatures of 
relativistic magnetic monopoles. An extensive collection of monopole references covering theoretical investigations 
and experimental searches prior to 1998 can be found in [3]. 

The electromagnetic interaction cross-section for relativistic monopoles is fairly well understood, but the hadronic 
cross-section is not. We calculate in some detail the reliable signatures resulting from the electromagnetic monopole- 
matter interaction. For the hadronic interactions of the monopole, we investigate various possibilities, and we present 
qualitative arguments for the resulting signatures. 

On theoretical grounds, the flux of monopoles that do not catalyze baryon decay is limited only by Parker's upper 
bound Fp ~ 10~^^/cm^/s/sr [4]. We do not consider monopoles which do catalyze baryon decay, as the fltix limits 
on these render them not observable. The Parker bound results from requiring that monopoles not short-circuit our 
Galactic magnetic fields faster than the Galactic dynamo can regenerate them. The Parker bound is several orders of 
magnitude above the observed highest-energy cosmic ray flux. Thus, existing cosmic ray detectors can meaningfully 
search for a monopole flux, and proposed vast-area cosmic ray detectors may improve the search sensitivity by many 
orders of magnitude. 

Because of the small inelasticity of monopole interactions, the monopole primary will continuously induce an air- 
shower. This is in contrast to nucleon and photon primaries which transfer nearly all of their energy in the shower 
initiation. Thus we expect the monopole shower to be readily distinguished from non-monopole initiated showers. 
However, we also investigate the possibility that the hadronic interaction of the monopole is sufficiently strong to 
produce air-showers with dE/dx comparable to that from nuclear primaries, in which case existing data would 
already imply a meaningful limit on the monopole flux. One may even speculate that monopoles with a large dE/dx 
have been observed, as the primaries producing the enigmatic showers above the GZK cutoff at 5 x 10^^ eV [5,6]. 

The outline of this paper is as follows: In section (II) we discuss the limits on the monopole mass and number density. 
We discuss the variety of monopoles resulting from possible symmetry breakings, and examine their natural kinetic 
energies resulting from acceleration in large-scale cosmic magnetic fields. In section (III) expressions for the energy 
loss of relativistic monopoles are presented. The electromagnetic processes of ionization, electron-pair production, 
bremsstrahlung, and the photonuclear interaction provide the most reliable expressions. The hadronic energy loss 
is also considered. In section (IV) we discuss signatures of relativistic monopoles arising from their electromagnetic 
energy loss processes. We develop a model for the electromagnetic particle shower induced by relativistic monopoles. 
Cherenkov signatures are also examined, including the coherent radio- Cherenkov signal resulting from the shower 
charge excess. We mention the possible detection of the coherent Cherenkov signal in km'^ ice experiments. It is 
shown that for certain masses Earth tomography with monopoles is practicable. The monopole's electromagnetic 
dE/dx in air appears to be too low to allow their detection via air fiuorescence. However, we argue that some 
speculative models may allow the geometric size of "baryonic-monopoles" to grow as a result of interactions, with the 
concomitant increase in their strong cross-section leading to significant air fluorescence. This possibility is explored 
in section (V) with a simple model for the growing strong cross-section. Finally, in section (VI) we give concluding 
remarks. 
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II. MONOPOLES AND MAGNETIC FIELDS 



A. Monopole Mass and Number Density 

As first shown by 't Hooft and Polyakov [7], the dynamical content and stabiUty of magnetic monopoles in particle 
physics models is regulated by the patterns of symmetry breaking in these models. A classification is given in [8] 
based on homotopy theory. For illustration, we consider the simplest case of a semisimple gauge group G (i.e., G has 
no f/(l) factors) that breaks to a subgroup G' that is nonsemisimple (i.e., has at least one U{1) factor which is part 
of a linear combination making up Uem{^))- At the symmetry-breaking temperature Tc {^) at which the order 
parameter {^) turns on, monopoles of mass M Tc/a appear, where a is the fine structure constant at scale Tc- The 
order parameter can be the VEV of some scalar field or some bi-fermionic condensate that breaks the symmetry. We 
assume that the value of (0) is at or above the electroweak (EW) scale 250 GeV so as to avoid violations of Standard 
Model (SM) physics. This assumption then bounds the monopole mass for our consideration to be ~ 40 TeV. We 
note that a similar lower mass bound results from (boldly) treating a classical monopole as a virtual quantum in 
radiative corrections to the SM [9]. 

The number density and therefore the flux of monopoles emerging from a phase transition are determined by the 
Kibble mechanism [1]. At the time of the phase transition, roughly one monopole or antimonopole is produced per 
correlated volume. The resulting monopole number density today is 

UM ~ 10-19 {n/10^^GeVf{lH/^cf cm-3, (2.1) 

where is the phase transition correlation length, bounded from above by the horizon size Ih at the time when 
the system relaxes to the true broken-symmetry vacuum. In a second order or weakly first order phase transition, 
the correlation length is comparable to the horizon size. In a strongly first order transition, the correlation length is 
considerably smaller than the horizon size. In section (II C) we will show that free monopoles with M ~ 10^^ GeV 
are accelerated to relativistic energies by the cosmic magnetic fields [10]. From Eq.(2.1) then, the general expression 
for the relativistic monopole flux may be written [6] 

Fm =cnM/47T ~ 2 X 10-" ( iQif^ev ) (f") ^m-" ^ec-^ sr-^ . (2.2) 

Phenomenologically, the monopole flux is constrained by cosmology and by astrophysics. Cosmology requires that 
the monopole energy density Om not be so large as to add observable curvature to the Universe. From eq. (2.1) comes 
an expression for the monopole mass density today relative to the closure value 

Qm ~ 0.1 {M/10^^GeV)\lH/^cf. (2.3) 

Monopoles less massive than - 10i3(^JZ^)3/4 GeV do not over-curve the Universe. Astrophysics requires the 
monopole flux to respect the Parker bound such that the magnetic field of our galaxy is sustained. Requiring that the 
Kibble flux in eq. (2.2) be less than the Parker limit Fp = 10-^^/cm^/sec/sr, one derives a combined mass bound [6] 

M ^10^\^c/lH)GeV . (2.4) 

This constraint is stronger than the curvature constraint by about two orders of magnitude for phase transitions 
without excessive latent heat. 
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The energy-density constraint for relativistic monopoles is of course stronger than that for non-relativistic 
monopoles of the same mass. Prom eqs. (2.3) and (2.2) we may write for the relativistic monopole closure density [6] 



where mpi = 1.2 x 10^^ GcV is the Plank mass. This shows that a Kibble monopole flux respecting the Parker limit 
cannot over-curve the Universe regardless of the nature of the monopole-creating phase transition (parameterized by 
^c/Ih), as long as (Em) ^ mpi. 

Minimal SU{5) breaking gives monopoles of mass ^ 10^^ GeV. However, within field theory there exist many 
possibilities to produce monopoles with mass below this scale. For example, there are other chiral SU{N) [11] and 
chiral 0{N) models [12,13] with lighter monopoles. There are also relatively light monopoles (M ~ 10* GeV) in 
the phenomenologically interesting SU{15) model [14]. Recently, a new field theoretic possibility has emerged based 
on conformal field theory (CFT). Although originally conceived in the context of type IIB strings and M- theory 
compactifications on higher-dimensional anti-de Sitter space [15], the rules for constructing the four-dimensional 
CFTs can be considered independently of any higher dimensional origins. These theories are naturally N = A 
supersymmetric, but depending on the compactification, can have reduced AT = or 1 supersymmetry and have 
sensible phenomenology with a product of SU{Ni) gauge groups. Three-family non-supersymmetric models have 
been constructed with GUT scale as low as a few TeV [16]. These models have light monopoles 100 TeV. Some of 
these models unify Uem{^) and SUc{^) [17], while others do not. Thus, monopoles in these models may or may not 
have hadronic interactions. 

In yet another alternative, extra dimensions a la Kaluza-Klein [18,19] decrease Mqut if these extra dimensions 
are not too tightly compactified. For example, with the introduction of two extra dimensions of millimeter size, 
the GUT scale is lowered to about 100 TeV due to the dramatic change from logarithmic to power law running in 
the RG-improved running coupling constant. This low-scale unification in turn may lead to magnetic monopoles of 
mass ~ 10'* TeV. A potential concern for monopoles in higher-dimensional theories is topological instability, i.e., 
monopoles in greater than four dimensions may unwind. However, monopoles made of SM fields confined to our 
three-dimensional brane would remain stable, In fact, there are even richer possibilities for monopoles in higher- 
dimensional theories. New defect solutions may exist in the larger dimensions, stabilized by topology. The intersection 
of these new solutions on the SM brane could then be monopoles, stabilized by the usual 3D topological argument. 

A relatively low GUT scale seems to be necessary to allow light mass monopoles. A low GUT scale will typically 
also enable fast proton decay, in violation of experimental lifetime bounds. However, exceptions abound. The SU{lb) 
model has no direct gauge induced proton decay. Proton decay in the MSSM Pati-Salam Model [11] is only induced 
through the hidden sector. Some of the 4D models listed above have accidental symmetries (continuous or discrete) 
which stabilize the proton, and the higher dimensional models have additional stabilizing possibilities. 

To summarize this section, there are theoretical possibilities for producing monopoles with mass ~ 10^^ GeV, 
while avoiding proton decay. In addition to their electromagnetic properties, these monopoles may have a strong 
interaction cross-section. In the context of the Kibble mechanism for monopole production, observational bounds on 
the Universe's curvature constrain the monopole mass to less than 10^^ GeV. More constraining is a comparison of 
the Kibble flux to the Parker limit which constrains the monopole mass to less than 10^^ GeV. 

We note that in higher dimensional cosmologies the Kibble flux given in eq. (2.2) may be altered. If the Kibble flux 




(2.5) 
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estimate is changed, then the straightforward Parker upper Umit Fm < 10~^^/cm^/sec/sr becomes the only rehable 
bound on the monopole flux. Thus, in the spirit of generality, we will let M be a free parameter and use the Kibble 
mechanism as a rough guide to Fm- We will, of course, require that Fm obey the Parker limit. We also will assume 
that proton decay is avoided in a way that does not restrict the parameter M. 



The fact that monopolcs arc topological defects endows them with a non-trivial internal structure. Monopoles are 
classified [1] by their topological winding in the group manifold. This topological classification is coarse for GUT 
monopoles, which require further classification according to their charges. Monopole charges are dual to the SM 
charges. For example, in a dual theory magnetic charge and chromomagnetic (or color-magnetic) charge replace 
electric charge and chromoelectric (or color-clcctric) charge. 

Fundamental monopoles can bind to form composite monopoles. For example, in an SU{b) GUT the fundamental 
minimally-charged monopoles are six-fold degenerate. For an appropriate Higgs potential there are four other types 
of stable bound states formed from the colored fundamental monopoles [20] . Remarkably, the spectrum of bound- 
state monopoles corresponds almost exactly to the particle spectrum of a standard model family, which has led to an 
attempt to construct a "dual standard model" out of monopoles [21]. 

We adopt the nomenclature "g-monopoles" for those monopoles with color-magnetic charge, and "Z-monopoles" 
for those with only the ordinary Uem{'^) magnetic charge. The possible confinement of q- monopoles via the formation 
of color-magnetic "strings" has been considered recently [22]. If such a confinement mechanism is realized, one 
expected result would be the formation of color-singlet "baryonic-monopoles." Thus we are led to consider the 
phenomenology of two broad classes of monopoles, the Z-monopoles with U{1) magnetic charge, and the baryonic- 
monopoles with hadronic interactions in addition to electromagnetic interactions. 

The baryonic-monopole structure is quite different from that of an Z-monopole, and as such we expect it to have 
a very difl[erent cross-section and cosmic ray shower profile. The internal dynamics of a baryonic-monopole would 
approximate that of an ordinary baryon in the QCD string model, but with g'-monopoles replacing quarks at the 
ends of strings. In particular, the string tensions of the chromomagnetic and chromoelectric strings are both of 
order A-qcd- However, the energetics of string-breaking for the two cases are very different. The chromoelectric 
string of standard QCD will stretch until energetics favors breaking and the production of a quark- ant iquark pair. 
The chromomagnetic string internal to the baryonic-monopole may stretch until energetics favors breaking and the 
production of a monopole- antimonopolc pair. From the large mass of the monopolc-antimonopolc pair, we infer that 
an enormous amount of string-stretching without breaking is possible for a baryonic-monopole. In sec. (V) we discuss 
the viability of baryonic-monopoles as candidate primaries for the super-GZK airshowers. 



The kinetic energy imparted to a magnetic monopole on traversing a magnetic field along a particular path is [6] 



B. Monopole Structure 



C. Monopole Acceleration 




path 



(2.6) 
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where 



g = e/2a = 3.3 x 10"^ esu (or 3.3 x IQ-Mynes/G) (2.7) 

is the magnetic charge according to the Dirac quantization condition, B is the magnetic field strength, ^ specifies the 
field's coherence length, and ^Jn is a factor to approximate the random-walk through the n domains of coherent fields 
traversed by the path. In Table 1 we indicate the cosmic magnetic fields [23] and their coherence lengths, inferred from 
observations of synchrotron radiation and Faraday rotation, and from modeling. The typical values of the monopole 
kinetic energies that result from these fields are also shown in the Table. A similar table has been constructed in ref. 
[24]. 

If the early universe dynamics that generated these fields and/or the present dynamics that maintain these fields 
were known, then new Parker like bounds could be placed upon the cxtragalactic monopole flux. With our limited 
knowledge at present, efforts in this direction (latter references in [4]) are highly model dependent. 

The largest energies arc seen to come from the magnetic fields having the longest coherence lengths. The strength 
of these sheet fields is inferred from simulations [25] and observations [26]. It is anticipated that in the near future 
more reliable inferences of the extra-galactic magnetic fields will become available [28,31]. This will allow a firmer 
prediction of monopole energies. 

We emphasize that a typical monopole which travels through the Universe, and has a mass below the energies 
indicated in Table 1, should be relativistic. Monopoles will gain and lose energy as they random-walk through the 
Universe, eventually producing a broad distribution of energies (with AEk/Ek ^ 1) centered roughly on ^/n times 
the typical energy for a single transit through a region of homogeneous magnetic field. Here, n is the number of 
coherent fields encountered in the random walk. For cxtragalactic sheets, which we expect to dominate the spectrum, 
this number can be roughly estimated to be of order n if(j'^/50Mpc 100. The resulting £^max is therefore 
estimated to be 10^^ eV. Hence, monopoles with mass below 10^^ GeV are expected to be relativistic. This is 
a fundamental result. The rest of this paper is devoted to the novel phenomenology of relativistic monopoles. We 
begin with a discussion of the interactions of monopoles with matter in the next section, and subsequently calculate 
monopole signatures in various detectors. 







C/Mpc 


gB^/eV 


Refs. 


normal galaxies 


3 to 10 


10-2 


(0.3 to 1) X 10^1 


[27] 


starburst galaxies 


10 to 50 


10--^ 


(1.7 to 8) X 10^" 


[28] 


AGN jets 


~ 100 


lO-"* to 10-2 


1.7 X (102" to 10^2) 


[29] 


galaxy clusters 


5 to 30 


10"* to 1 


3 X 10^* to 5 X 102^ 


[30] 


Extragal. sheets 


0.1 to 1.0 


1 to 30 


1.7 X 10^^ to 5 X 1023 


[25] 



TABLE I. Estimated magnetic field strength and coherence length for some astrophysical environments, and the associated 
monopole energies for a single transit through the regions. 



6 



III. RELATIVISTIC MONOPOLE ENERGY LOSS IN MATTER 



GUT monopolcs arc formed in non-trivial representations of SUc{'^) x SUl{'2) x C/y(l) and couple to the standard 
model gauge fields. An accurate description of the monopole stopping power must include all of the standard model 
interactions or the relevant subset of those interactions for the type of monopole considered. In this paper we choose 
to ignore the weak interaction throughout, which is suppressed in amplitude by factors ~ M^^. The electromagnetic 
and strong interactions will be accounted for in this section. 

The strong interaction of a monopole is difficult to assess. Color confinement ensures that all observable monopoles 
are color singlet objects, and so have no classical long range color-magnetic field. Nevertheless, we expect Z-monopoles 
and baryonic-monopoles to have very different hadronic interactions. Although Z-monopoles are fundamental and 
lack a color-magnetic charge, the unbroken symmetry in their core ensures that gluon and light quark fields will 
leak out from the center to the confinement distance ^qcd ^ 1 fm [32]. In this way, 1- monopoles will exhibit 
some hadronic cross-section in small impact-parameter scattering. This is probably ignorable. On the other hand, 
baryonic-monopoles are intrinsically hadronic in all partial waves. 

We will resume the discussion of the monopole's strong interaction with hadronic matter after first discussing in 
some detail their better-understood electromagnetic interactions. The electromagnetic interaction of the monopole 
may dominate the hadronic interaction because the electromagnetic coupling of the monopole is large [33], 

aM = ^^M, (3.1) 

and mediated by a long-range field. At large distances and high velocities, the magnetic monopole with Dirac charge 
mimics the electromagnetic interaction of a heavy ion of charge Z ^ \JoimIol ^ ^ — 68. We will follow others and 
treat the electromagnetic interaction of the monopole with matter semi-classically, i.e., viewing the monopole as a 
classical source of radiation, while treating the matter-radiation interaction quantum-mechanically. In this way, the 
large electromagnetic coupling of the monopole is isolated in the classical field, and the matter-radiation interaction 
can be calculated perturbatively [34]. 

The value of au given here, and conservatively used throughout this paper, is the minimal value due to Dirac. 
Other monopole solutions have charges which are integer multiples of the Dirac charge, and therefore ckm's which 
are larger by a factor of . n = 1, 2, 3, . . .. For example, the Schwinger solution and the original t'Hooft-Polyakov 
monopole have twice the Dirac charge, and an q;m that is four times minimal. 

A. Electromagnetic Interactions 

We consider here the energy losses of the monopole due to the four electromagnetic processes: ionization resulting 
from collisions with atomic electrons, e+e~ pair production, bremsstrahlung, and the photonuclear interaction. It 
will be the pair production and photonuclear interactions, that dominate the energy loss for very large 7 = E/M and 
will be mostly responsible for the growth of a particle shower. The monopole-matter electromagnetic interaction for 
7 < 100 is well reported [35,36] for atomic excitations and ionization losses in the absorber, including the density- 
dependent suppression effect arising from charge-polarization of the medium [37]. In the literature these two effects 
are collectively referred to as "collisional" energy losses and we follow that nomenclature. 
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It is convenient to express the energy loss per "length" in units of column density (g/cm^) which we use throughout 
this work. The conversion from true length increment to column density increment dx is dL = dx/pM = 
(^)(^) X, the latter form arising from writing the nuclear matter-density p as tun nN, and identifying the nuclcon 
mass niN and density un with the first and second fractions, respectively. In conventional notation, Na is Avogadro's 
number, Z and A arc the mean nuclear charge and number, rig is the electron density, and g is the gram unit. We 
absorb the mass into A, giving A the units of grams per mol. 

One approach to electromagnetic interactions is to replace the monopolc current with an equivalent photon flux, 
dNy/duj = 2Q:M/7rti^ lii(TO7/a;), with m being the target mass and lu the photon energy in the rest frame. The 
equivalent photon approximation (EPA) is strictly valid for limited perpendicular momentum transfer q± <^ to, and 
limited energy transfer uj <C mj [38]. Consequently, we prefer to use the more exact monopole-target scattering 
formula. Nevertheless, the EPA presents two illuminating features, the l/u bremsstrahlung-type of photon spectrum, 
and the large normalization 2aM7r = l/27rQ! w 22. The latter tells us that electromagnetic cross-sections are large, 
while the former, in conjunction with u) <C rwy, shows that energy transfers are soft. 



1. Collisional Energy Loss 

The magnetic monopole stopping power formula calculated by Ahlen [36] includes ionization of the absorber at 
small impact parameters and atomic excitation at large impact parameters. For highly relativistic monopoles we may 
ignore the various correction terms and simply describe ionization with 

dEcou _ 4:nZ(\(\MNA 
dx ArUe 

me is the electron mass, / is the mean charge and ionization energy of the material, and 6 parameterizes the density 
effect [39]. 



In 



(3.2) 



2. Pair Production 

A thorough study of pair production for muons [40] is adapted here for the case of a monopole. To accurately 
calculate shower development resulting from pairproduction the inelasticity must be understood. For a relativistic 
monopole primary the energy loss to pair production is 

where 



TOfi 



Ex- (3-4) 



Xi= f''"'^^ driFi{-f,ri), (3.5) 



the inelasticity, or fraction of monopole energy transferred in the interaction, is ?7 = AE/E, and the five functions 
Fi{'y,r]) (given in [40]) are labeled by the five distinct regions of inelasticity and nuclear-charge screening: 
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1) slow pairs, no screening, < r] < 



2) slow pairs, total screening, ^ < rj < ^ 



jM ^ '/ ^ M 

7M ^ '/ ^ M 

3) fast pairs, no screening, < r] < (3.6) 

4) fast pairs, no screening, <r} <1 

me 

5) fast pairs, total screening, < r/ < 1. 




FIG. 1. Electron-pair production of an M = 100 TeV relativistic monopole in air is described in terms of the functions Xi 
derived in [40] . The Xi 's relate to different regions of inelasticity space and either total or no screening. 



In fig. (1) wc plot the factors versus 7 to examine the relative strengths of the Xi's- To a good approximation, 

pair production is proportional to the sum of just X2 and X3. The expression xi is unphysical over the range of 7 
plotted while Xi and Xb have been multiplied by 10^. For M ^ me, and keeping mass corrections with the largest 
powers of ^ In" (^^) ; the no screening limit is 



pair 



dx 



197r a^amZ'^NK 
9 Anie 



(l-Bi)lnQ)-B2' 



and the total screening limit is 



dx 



197r oF'aMZ'-Np^ 
9 Ame 



(l-Si)ln 



189 
ZV3 



In 2 - S3 



(3.7) 



(3.8) 



where 



R — 48 m. i„2 ( m\ 
R _ 11 16 me 1^3 / M 

^2 - -g- - 19^ irin I — 



D 1 _|_ 16 me 1^^ 

^3 - 38 + WW'U^ 



(3.9) 
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For E ^ M ^ rUe Bs here, the contributions of the three Bi to eqs. (3.7) and (3.8) are negligible. Examination 
of fig. (1) shows that for 7 ^ 10^, a total screening process dominates, meaning that eq. (3.8) describes the pair 
production energy loss. 



3. Monopole Bremsstrahlung 



Bremsstrahlung radiation by a rclativistic particle in collision is inversely proportional to its mass, so we expect 
bremsstrahlung to be negligible for massive monopoles. The approximate stopping power due to brchmsstrahlung 
losses has been calculated [41] for electric charges in the negligible nuclear-screening limit [42]. For monopoles it is 



4- Monopole Photonuclear Interaction 



Virtual photon exchange between a monopole and a nucleus is described by the photonuclear cross section. This 
process is MN MX, where M is a monopole, N a nucleus, and X are final state hadrons. The energy loss of 
Icptons via the photonuclear interaction has been re-cvaluated recently [43] in light of the HERA results for real and 
virtual photon nucleon scattering. A full range of momentum transfer is accounted for from real (Q 0) photon 
exchange to highly virtual photon-nucleon interactions. This analysis is adapted to monopoles with the replacements 
rUfj, M and e e/2a where 171^, M are the muon and monopole masses, respectively. For 7 > 10^ the photonuclear 
energy loss dominates pair production and the calculations, based on [43], give roughly dE/dx oc 7^-^®. 



B. Total Electromagnetic Losses 



We collect the electromagnetic energy loss processes together and plot them in fig. (2) for M = 100 TeV monopoles. 
The dominant interactions are atomic collisions for 7 < 10^, pair production for 10^ < 7 < 10^, and photonuclear for 
7 > 10^ 
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E 



M = 00 TeV 
Air, 1=1.1 




o 



Boost Factor, 7 



FIG. 2. The electromagnetic energy loss from collisional, bremsstrahlung, electron-pair production, and the photonuclear 
interaction of a 100 TeV relativistic monopole in air. Collisional, pair production, and the photonuclear interaction are roughly 
independent of the monopole mass wherecis bremsstrahlung is ex . The units of energy loss are given in TeV per atmosphere, 
where 1 atm = 1033 g/cm^. 



GUT monopoles typically contain internal color fields. The l-monopole's color field is soft, extending to a distance 
^ 1 fm (where A = Kqcd)- Its internal color is similar to that of a normal hadron, however there are no 
valence quarks internally, and so a better analogy to the color structure of an Z-monopole is a glueball. Thus, we 
expect the hadronic cross-section of Z-monopoles to be of typical hadronic size, with an inelasticity r\ ~ A/M ~ 
10~^° (10^'^GeV/M). For Z-monopoles an approximate hadronic stopping power would then be 



where the mean-free-path between hadronic interactions is A = (iVAUhad)"^, for a typical hadronic cross-section 
o'had- The magnitude of expression (3.11) appears to be comparable to the energy loss to electron-pairs discussed 
above, so the hadronic interaction for Z-monopoles may make a significant contribution to the total energy loss. 
Lacking reliable knowledge of the Z-monopole's hadronic interactions, we will neglect this possible contribution in the 
following signature calculations. In passing we mention that: 1) The inclusion of eq. (3.11) in the longitudinal shower 
model (sec. IV) would feed the electromagnetic shower. For each hadronic interaction length in the development of a 
subshower, about 30% of the hadronic energy is transferred to the electromagnetic component; thus, after a subshower 
undergoes a few hadronic interactions, 0(1) of the energy will be in the electromagnetic shower. 2) The range of mass 



C. Hadronic Interactions 




(3.11) 
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for which monopole tomography (to be discussed in in sec. (IV D)) is viable will be shifted toward higher mass with 
the inclusion of eq. (3.11). 

The baryonic-monopoles are quite different. In section (V) we develop a model for the cnca'gy loss of baryonic- 
monopoles. The QCD string model, where the total cross-section grows with string length /, describes baryonic- 
monopolc hadronic interactions, with the significant caveat that the confining strings in the monopole are readily 
stretched but not easily broken. The energy loss is estimated to be 

(^) - -TTT - -7 A cr(a;) ~ -7 l{x) , (3.12) 



dx X{x) \gram/ \gram 

where the string cross-section a{x) = l{x)A~^ is explicitly a function of column depth x 



IV. MONOPOLE ELECTROMAGNETIC SIGNATURES 



Signature events for Z-monopoles are derived below with a specific emphasis on 1) the general shower development, 
2) Chcrcnkov signatures, 3) the nitrogen fluorescence signature, and 4) Earth tomography. The general shower 
characteristics are developed first as the other signatures arc derivable from that model. For the remainder of this 
section we only consider /-monopoles, which we will simply refer to as "monopoles" . 



A. Monopole Shower Development 



Monopoles will be highly penetrating primaries. On average, there will be a quasi-steady cloud of secondary particles 
continuously regenerated along the monopole trajectory. Thus, we will call this type of shower "monopole-induced," 
in contrast to conventional particle-initiated showers. 



1. Monopole-induced Subshowers 



Given a fast monopole passing through matter, the various electromagnetic processes discussed in sec. (Ill) can 
inject energetic photons, electrons, and positrons into the absorbing medium. If the energy of these injected secondary 
particles is sufficient, they may initiate a particle cascade. In the simple model we consider, originally developed by 
Heitler [44] and reviewed in [45,46], the photon pair production length is set equal to the electron (or positron) radiation 
length. Consider a photon primary of initial energy Eq which travels a mean distance R through the absorber before 
pair producing [47] . The two particles in the produced pair are assumed to share equally the initial photon's energy. 
After traveling another distance R the electron and positron each radiate a bremsstrahlung photon where the produced 
photon takes half the parent energy. The particle number has doubled again and the energy per particle is halved 
again. The shower continues to develop in this geometric fashion until the energy per particle drops to the critical 
energy, Ec, below which ionization of the absorber dominates the energy losses due to particle-production processes. 
For typical materials, Ec ~ 100 MeV. 

From the relation for energy per particle, 

E{x) = Eoe-k, (4.1) 
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it is apparent that R is related to the radiation length by = lii2. The number of particles in the shower at a 
given depth is approximately 

N{x) = eft . (4.2) 
The maximum number produced is A''inax = ^ , from which we infer the depth of the shower maximum to be 

X^ax = eeln(^^) . (4.3) 

After reaching the shower maximum at Xmax, the shower size decreases with column depth as a rclativlcy mild 
exponential. The value of in air is quoted as 34g/cm^ [46], while the value in ice is quoted as 36g/cm^ [50]. The 
slight dependence on composition reflects the small differences in mean charge of the nuclei. The attenuation length 
for the post-Xmax decay is approximately 200 [51], or about 6^e- 



2. Longitudinal Shower Profile 

It is straightforward to qualitatively describe the electromagnetic shower from pair production. The subshowers 
begin their existence as an e+e~ pair, produced nearly at rest in the cms. The subshower energy is therefore ~ 2me 
boosted to the lab frame. Since M is much greater than the nucleon mass, the cms frame is nearly coincident with 
the monopole rest frame, and we have Eq 2'yme for the subshower energy. The electromagnetic process by which 
the subshower repeatedly cascades, and the process by which the monopole initiates the subshower differ only by the 
attachment of the monopole to a photon line; i.e., the monopole interaction mean free path A is a factor of ^ ^ 16 
smaller than the subshower radiation length ^e- Each subshower attains its A^max and retains this particle number for 
roughly an attenuation length 6 (as presented above). Thus, the number of subshowers contributing to the steady- 
state shower size at any given moment is 6^e/A ~ 100. With ^ 100 subshowers contributing to the total shower, we 
find N ^ WOEq/Ec ^ 7 for the total particle number of the steady-state cloud accompanying the monopole. It is a 
pleasant coincidence that the proportionality in A' ~ 7 is of order unity. 

The electromagnetic shower produced by the photonuclear process is harder to estimate. This process dominates 
the total energy loss for 7 > 10® and will dominate the fluorescence yield over that range, but it's contribution to the 
electromagnetic shower is indirect. In an Appendix we present a quantitative model for the electromagnetic shower 
from pair production and discuss the contribution from the photonuclear interaction. In fig. (8) shown later, we plot 
the shower size A'^max versus 7 based on the detailed shower model. It is reassuring to see the qualitative agreement 
between the detailed calculation and the approximate relation A^^ax ~ 7 with its pleasing 0(1) proportionality. 



B. Monopole Cherenkov Signatures 



When a charge travels through a medium with index of refraction n, at a velocity (3 > ^, Cherenkov radiation is 
emitted [52] . The Frank- Tamm formula gives the total power emitted per unit frequency v and per unit length I by 
a charge Ze, 



dudl 



1 - 



(4.4) 
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The maximal emission of the Cherenkov light occurs at an angle ^max = arccos(l/n/3) where 6 is measured from the 
radiating particle's direction. 

1. Direct Monopole Cherenkov Emission 

The interaction of a magnetic charge with bulk matter requires the replacement of factors of e with the Maxwell 
dual factors /i. But /i and e are related by the index of refraction. The replacement the electric charge formulae, 
eq. (4.4) with Z = I, adequate for magnetic monopoles is a ^ n'^/Aa, and leads to an enhancement factor of 4700 
for monopoles interacting in vacuum and 8300 for monopole interactions in water. Cherenkov light from an electric 
charge source is linearly polarized in the plane containing the path of the source and the direction of observation. 
However, the light polarization from a magnetic charge will be rotated 90 degrees from that of an electric charge 
which, in principle, offers a unique Cherenkov signature for monopoles [53]. 

The direct Cherenkov signature of a monopole leads to the best experimental limits at present on the flux of 
relativistic monopoles. Prom the absence of a Cherenkov signal in deep water PMTs the Baikal Collaboration reports 
the limit Fm ^ 6 x 10~^®/cm^/s/sr [54], while the absent signal in deep ice translates into the AMANDA limit 
Fm ~ 1.6 X IQ-i'^/cmVs/sr [55]. There is a comparable limit of ~ 4 x IQ-^Vc^niVs/sr obtained by the MACRO 
collaboration using scintillators and streamer tubes [56]. Note that these limits are already slightly more restrictive 
than the Parker limit. This work shows that the direct Cherenkov signal in the optical should be swamped by the 
electromagnetic shower for monopoles with 7 > 10^ (see fig. (6)). Such a signal is not seen by AMANDA [57] and so 
the AMANDA monopole flux limit remains the same as that derived from the non-observation of a direct Cherenkov 
signal. In sec. (IV D) we calculate that relativistic monopoles can pass through a large portion of the earth for M ~ 3 
PeV. Given this result, any experimental limit for upgoing relativistic monopoles should only apply roughly in the 
mass range M ~ 3 PeV. 

2. Coherent Radio-Wavelength Cherenkov Emission 

In addition to the Cherenkov radiation from the bare monopole charge there is a contribution from the relativistic 
e+-e~ cascade comprising the shower. Therefore, an estimation of the total power radiated in Cherenkov light requires 
the model developed above for the size of the monopole induced shower. 

The approximate lateral width of a monopole-induced cascade is given by the Molicrc radius (7.15) which means the 
excess electric charge [58] is confined roughly within a distance Rm- Cherenkov light of wavelength A ^ _Rm emitted 
by the monopole-induced shower will be radiated coherently and be strengthened by the factor in eq. (4.4). 
Simulations of electromagnetic showers in ice [50] show that the cascade contains an electric charge excess of about 
20% the total shower size. This is plotted for ice in figure (8). 

The proposed Radio Ice Cherenkov Experiment (RICE) may offer the best chance for the detection of monopoles 
in a Ikm'' scale detector resulting from the transparency of ice in the radio. [59] A thorough study, on a par with 
that done for high energy neutrinos [60], has not been undertaken for monopoles in RICE. In lieu of such a study we 
remark that the most noticeable signature should be that the monopole is highly penetrating and will traverse the full 
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detector size without an appreciable loss of kinetic energy. Thus, the monopole signature is easily distinguished from 
a neutrino event, which is a localized shower that produces a Cherenkov cone detectable only by a limited number of 
antennae lying within the cone. Timing can be used to reconstruct the monopole path. 

The analysis of neutrino energy thresholds at RICE can be used to estimate an energy threshold for monopole 
events. Figure (13) in [60] can be interpreted for monopole detection as showing that an effective volume of Ikm'' 
per antenna is reached if dE/dx ^ dEth/dx w 3 x 10^2 cVcm^/g « 3 X 10-'' TcV/atm. This translates into a threshold 
boost factor, read off from figure (2), of 7 ^ 7th ~ 10^. For a RICE monopole search to have an effective volume 
~ 1km the monopole boost is restricted to 7 ^ 10 . 

The monopole flux limits which RICE could set are much stronger than those currently available (Baikal, AMANDA, 
or MACRO [54-56,61]). A year of observation in the km'^ RICE detector without a candidate monopole event would 
set a flux limit at least as good as 10~^^ cm~^sec~^sr~^, three orders of magnitude below the Parker bound and two 
orders of magnitude better than the current bounds. 



C. Monopole Fluorescence Signatures 



A sensitive probe of high energy air shower development is the nitrogen fluorescence emitted by the atmosphere. 
The Fly's Eye and HiRes experiments have fruitfully pioneered this technique. The yield in nitrogen fluorescence is 

0.5% times the total energy in ionization of the atmosphere. Since 0{1) of the energy loss by all processes ends up 
in ionization of the atmosphere once the subshowers have ranged out, we need to evaluate the total stopping power 
to estimate the fluorescence yield. The total energy loss is the sum of the various processes summarized above: 

dx dx dx dx dx dx 

For Z-monopoles the electromagnetic energy loss dominates, and in a thin absorber like the atmosphere, Z-monopole 
energy deposition is typically below threshold [62]. For light Z-monopoles at extreme 7, photonuclear losses dominate 
and 

^ « - 10« ^ r A V f lilU^) (4.6) 
dx cosOz \7.2J \ A JKIO^J atm ^ ' 

On the other hand, baryonic-monopoles [22] interact such that dominates the total energy loss. For baryonic- 

monopoles a large energy transfer to the atmosphere is natural. In section (V) a model for the baryonic-monopole 
interaction is developed. 



D. Earth Tomography with Relativistic Monopoles 



Direct knowledge about the composition and density of the Earth's interior is lacking. Analysis of the seismic data 
is currently the best source of information about the Earth's internal properties [64,65]. However, another potential 
probe would be the study of highly penetrating particles interacting with the Earth's interior. With such a means it 
may be possible to directly measure the density profile of the Earth's interior [66] . Here we show that over a range of 
initial kinetic energies monopoles can pass through the Earth's interior and emerge with relativistic velocities. The 
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results of a numerical calculation, making use of the full energy loss expressions derived earlier and a simple model 

approximating the internal composition of the Earth, are described below. Our calculation shows that relativistic 
monopoles in the 3 ^ 10 PeV (1 PeV = 10^^ eV) mass range are ideal for Earth tomography, and that monopoles 
may pass through the Earth to initiate up-going cosmic ray events. 



1. Two Shell Earth Model 



Dziewonski and Anderson [67] have developed the preliminary Earth reference model (PREM) which is the standard 
Earth model in use today. The PREM model consists of eight concentric shells of varying density and composition. It 
is sufficient for the purpose of demonstrating monopole tomography to simplify this model to two shells, the mantle 
and the core [68]. Both shells are taken to be spherically symmetric. The core has a radius i?core = 3.486 x 10^ m and 
a mean mass density pcore = H-S^^- The mantle extends from the core out to the earth's surface at = 6.371 x 10^ 
m and has a mean mass density Pmantic = 4.0^^. We take both shells to be of uniform composition. The chemical 
composition of the mantle (in mass) is approximated by Si02 (45.0%), AI2O3 (3.2%), FeO (15.7%), MgO (32.7%), and 
CaO (3.4%), and that of the core is approximated by Fe (96.0%) and Ni (4.0%) [69]. From these data we calculate the 
chemical composition (in percentage of molecular type) for the mantle as Si02 (40.0%), AI2O3 (1.67%), FeO (11.6%), 
MgO (43.6%), and CaO (3.23%), and for the core as Fe (96.2%) and Ni (3.8%). 

The electromagnetic energy loss processes for relativistic monopoles were analyzed in sec. (HI). The energy loss 
formulae contain various parameters which need to be specified for the core and mantle of the Earth. The density 
effect includes the following dimensionful parameters: 
the plasma frequency [36] 

[ 66.2 eV = 3.36 x lO^^cm"! (core) 
I 40.4 eV = 2.05 x lO^^cm"! (mantle) 

and the mean ionization potential [36] 

' 285 eV (core) 
172 eV (mantle) 

and various numbers in the Sternheimer and Peierls parameterization [39] : 



(4.8) 



, / \ -3.92 (core) 

C=-2\n(—)-lc,{ ^ ' ), 4.9 

\<^pJ -3.90 (mantle) 



(7 + 4.606X0 _ I 0.137 (core) 
(ATi - Xof ~ j 0.136 (mantle) 



(4.10) 



the latter obtained from Xq = 0.2 and Xi = 3.0 for both the core and the mantle. The functional form [39] for 5(7) 
in the core is 



7 < 1.87 

jcore(^-) ~ ^ _3.92 + 111(72 - 1) + 0.0297 [13.82 - ln{-f^ - 1)]^ 1.87 < 7 < 1000 
^ -3.92 + ln(72 - 1) 1000 < 7 



} , (4.11) 
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and in the mantle is 

{0 7 < 1.87 

-3.90 + ln(72 - 1) + 0.0295 [13.82 - ln(72 - 1)] ^ 1.87 < 7 < 1000 } ■ (4.12) 

-3.90 + ln(72 - 1) 1000 < 7 

Lastly, the mean nuclear charges and mean atomic weights, calculated from the chemical compositions given above, 
are 

, 26.1 (core) , 

}, (4.13) 

10.8 (mantle) 

,56.1 grams/mol (core) , 
A~ { ^ ' } ■ (4-14) 

21.7 grams/mol (mantle) 



2. Numerical Results 

We integrate ^^-^^^ = jj — to find the final boost factor for the monopole after traversing the earth. The input 
parameters are the monopole mass and initial boost factor. The results are plotted in figure (4). The figure shows 
several curves with different initial kinetic energies, for a fixed monopole mass. 




FIG. 3. A schematic representation of earth tomography with a highly penetrating particle. Monopole energy losses differ 
in passing through the core or mantle and so affect the energy spectrum of upcoming monopoles as a function of zenith angle 
Oz ■ For a zenith angle 6^'^'^ ~ 147° , the upcoming particle grazes the edge of the core. 
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The initial energy spectrum can be determined from the < f data. These monopoles have not passed through 
the earth and therefore retain the original kinetic energy attained from the large-scale magnetic fields. A comparison 
of these data with the 0z > f data could make a determination of the earth's interior possible. In figure 4, the 
discontinuity due to the mantle-core boundary is clear [70]. The zenith angle for a trajectory tangent to the core 
boundary is 



Rc 

Ra 



'^ore = ^ _ arcsin ~ 147°. (4.15) 



Assuming that 0(100) data points would be sufficient to delineate the attenuation curves of figure 4 we can estimate 
the detector exposure needed in the 2n steradians projected through the Earth: 

1 / N \ /10-iVcm2secsr\ /lkm\^ 
t = lyear - — • (4-16) 



The monopole flux Fm ~ 10~^^/cm^ secsr is roughly the experimental upper bound of [54-56]. 





FIG. 4. M = 3 X 10^^ and M = 10^^ eV monopoles passing through the earth and emerging with a diminished boost factor 
for increased zenith angle. The range of initial kinetic energies presented here is 3 x 10^° eV < Em < 10^^ eV. 

V. BARYONIC-MONOPOLE AIR SHOWERS AND SUPER-GZK EVENTS 

The observation of air-showers above the GZK cutoff at Egzk ^ 5 x 10^® eV presents a puzzle. These events cannot 
be due to nucleons or photons propagating from sources located at cosmic distances (see [6] for a recent discussion of 
the puzzling nature of these events). The natural acceleration of monopoles to energies above the GZK cutoff, and 
the allowed abundance of a monopole flux at the observed super-GZK event rate of 2.0 ± 0.5 x 10~^°/cm^/s/sr above 
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10^" eV [71] (five orders of magnitude below the Parker limit), motivates us to ask whether monopoles may contribute 
to the observed super-GZK events. As a proof of principle, here we present a simple model of a baryonic-monopole 
interaction in air which produces a shower similar to that arising from a proton primary. (We note that an unrelated 
mechanism for mimicking a high-energy proton-initiated air shower with a monopole-nucleon bound state has been 
advanced in [24].) To mimic a proton-induced shower the monopole must transfer nearly all of its energy to the shower 
in a very small distance. The large inertia of a massive monopole makes this impossible if the strong cross-section is 
typical, ~ 100 mb [72]. The cross-section we seek needs to be much larger. 




(a) (b) (c) (d) 

FIG. 5. A schematic representation of a baryonic-monopole interacting with the atmosphere that depicts the effect of 
string-nucleon interactions. Fig. (a) shows the baryonic-monopole in its unstretched state, before hitting the atmosphere. 
After the first string-nucleon interaction (fig. (b)) the string stretches to length ~ Eeti A~^, and after n interactions (fig. (c)) 
it stretches to ~ n Eetr A^^ . Fig. (d) shows string stretching from g-monopole recoil. 

We model our arguments on [22] where three g— monopoles are confined by Z3 strings of color-magnetic flux to 
form a color-singlet baryonic-monopole as in figure (5a). All scales in the ground state baryonic-monopole are set by 
the QCD scale A. In particular, the string tension has the usual QCD-strength /i ~ A^. We further assume that: 

(1) Before hitting the atmosphere, the baryonic-monopole's cross-section is roughly hadronic, ctq ^ A^^. 

(2) Each interaction between the baryonic-monopole and an air nucleus transfers on average energy i?str (in the 
monopole rest-frame) into the chromomagnetic string system of the baryonic-monopole. 

(3) Since the chromomagnetic strings can only be broken with the formation of a heavy (7— monopole-antimonopole 
pair, the internal energy does not break the strings but rather stretches them by amount 6L = Est-c/ ^ ~ (i?str/A) fm; 
the string system grows linearly with each interaction, until other emission processes of the stretched (and presumably 
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vibrating) string dominate [73]. 

(4) The cross-section for the interaction of the baryonic-monopole with a nucleus is geometric, a ^ LA~^, where 
L is the (growing) string length; this assumption says that string-nucleon interactions dominate over g-nucleon 
interactions. 

(5) The energy emission rate from the excited string system (in the monopole frame) is proportional to some small 
power p of the (growing) string length (and therefore, to the same power of the number of hits and (growing) cross- 
section). 

Together, the assumptions imply that the baryonic-monopole cross-section grows from typical hadronic at first 
interaction, ctq ~ A~^, to 

an = LK~^=n^^ (5.1) 

after n interactions. A scematic drawing of the baryonic-monopole is given in figure (5) for (a) the ground state, 
(b) after the first interaction, and (c) after the n*^ interaction. This growth (ex L) continues with each interaction, 
until finally the rising rate of particle emission E'emit (oc U') equals the rate of energy acquisition due to stretching. 
If iJstr ^ A, then already after the Rrst interaction the cross-section is sufficiently large to shrink the subsequent 
interaction length to a small fraction of the depth of first interaction. This allows 0{1) of the incident monopole 
energy to be transferred to the air nuclei over a short distance, just as in a hadron-initiated shower. 

We can be more quantitative. The rate of internal energy acquisition eflfecting string-stretching exceeds the rising 
rate of particle emission as long as 

^>4mH, (5.2) 

where At is the time between interactions in the monopole rest frame. Explicitly, we have 

(Af)-i = 7 (Atiab)-' = 7t^ = c-ipN<J , (5.3) 

Alab 

where pjv is the nucleon density of air in the lab frame. It is important to the model that the monopole first has a 
growing phase, and then before it becomes non-relativistic, the emission phase. This requires that inequality (5.2) 
hold for the initial interactions. Such will be the case when p > 1, i.e. the emission rate from the string grows 
nonlinearly with the string length. We comment on this below. Inputting eqs. (5.1) and (5.3) into (5.2), we find that 
the baryonic-monopole remains in its growth phase until the critical n^^ interaction 



A2 

-C'str 



(5.4) 



A k{p) 

Here, k{p) is the proportionality constant in E'emit = k{p)U'. The "size" of the monopole at the conclusion of its 
growth phase is an^^, = ristr-Estr/A^. (For our order of magnitude estimates, we benignly neglect the recoil of the 
monopole and the related decrease of 7 with each interaction.) The distance between the first interaction and the 
onset of the emission phase is 



"str 

AX - ^ A„ , (5.5) 

n=l 



which, according to A„ = AAq / Egtrii, is 



20 



Ao and A„ arc the mfp's of the monopole in the initial ground state and after n interactions, respectively, and the 
density of nucleons is taken here to be a constant for simplicity of illustration. The sum is a finite series, and for 
ristr ^ 1 is very nearly equal to In(nstr) [75], giving 

AXstr ~ Ao In(nstr) (5.7) 

^str 

as the estimate for the total distance traveled by the baryonic monopole in its growth phase. As derived here, AX 
and A arc lengths in the monopole rcst-framc. However, our cqs. (5.5) (5.7) are homogeneous in these variables, and 
so also apply in the Earth frame without change. In terms of column density of air, this distance is just 

Aa;str = ^ AX ~ 100 ( ^) ln(n«tr) g/cm' , (5.8) 

where the numerical value is evaluated in the Earth frame. Turning to the emission phase, equilibrium between 
absorption and emission fixes the average energy emitted per interaction to be ~ E^str in the monopole frame. An 
energy M in the monopole frame, isotropically emitted, boosts to 7M = E m the lab frame. Therefore, the number 
of interations required to absorb and re-emit the energy observed as an air-shower is 

"emit ~ t;— • (5.9) 

-C'str 

Numbers are Lorentz invariants. Using Xn^^^ = Ai/rigtr = ^nd eq. (5.7), the distance traveled in the emission 

phase is then 

MA Ao 

This too can be quite short compared to the mfp Aq of the ground-state monopole. 

The two constraints on the model, that the length of the stretching phase AXgtr and the length of the emission 
phase AXemit be short compared to the mfp Ao for first interaction, can be written with the help of eqs. (5.7) and 
(5.10) as: 

^'^«n,t.«e^-/^, (5.11) 



AXemit ~ '^emit A„^^^ — 2 Z. " (5.10) 



p2 

-'^str 



with ristr given by cq. (5.4). For a monopole mass of 100 TcV (10^" GcV), the left and right terms of the inequahty 
arc correctly ordered for i?str/A > 10 (20). Thus, for iJgtr 3> A, a broad range of rigtr values will satisfy cq. (5.11). 
Furthermore, all values of rigtr are available by suitably choosing the free parameter Thus, the model is natural in 
that it does not require fine-tuning of parameters. We have succeeded in constructing a stretchable chromomagnetic- 
string model for the baryonic-monopole which provides an existence proof that a very massive monopole, despite its 
inertia, may nevertheless transfer 0(1) of its relativistic energy to an air shower over a very short distance. Such a 
monopole mimics the signature of a primary nucleon or nucleus. Of course, the last gasp of the monopole will not 
resemble a nucleon, for the monopole will eventually become non-relativistic and shed its internal energy more and 
more isotropically. If this "end game" energy is a small fraction of the air-shower energy, it will be difficult to observe. 

Let us discuss assumptions (2), (3), and (5). Any assumption that does not violate known physics is valid until 
disproven by Nature. However, it is motivational to look for plausibility arguments. Assumption (2) states that 
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energy Estr per interaction stretches the string. An analogy to p — p scattering may be illuminating here. The 
baryonic monopole contains a QCD energy of roughly A. Thus, the stong-interaction scattering dynamics bear some 
resemblance to p~p scattering at s ~ jAtun ~ (TeV)^. p — p scattering data at 2 TeV cms energy are available from 
Fcrmilab's p — p Tevatron collider. A mainly diffractivc QCD interaction stretches particle production throughout 
the rapidity plateau, emitting ^ 10^ particles with typically ^ few GcV energies. When a string is struck by 
an air nucleus, the diffractivc QCD dynamics may bo similar, stretching the colored strings across the rapidity 
plateau. However, by assumption (3), the chromomagnctic string cannot break, and so a significant fraction of the 
0(100 GeV) energy remains in the stretched string, with the remaining energy materializing as light hadrons. The 
value Estr ~ 0(100 GeV), implying _Bstr/A ^ 10'^, appears more than ample (by an order of magnitude or more) to 
sustain the model we have outlined. We also note that whereas the color-electric charge interactions are asymptotically 
free, the strength of the color-magnetic interactions increases with energy to preserve the Dirac quantization condition. 
This may further argue in favor of the non-breaking, coherent behavior of the chromomagnetic string. 

Assumption (3) states that energy emission by the excited string is negligible in the early growth phase of the 
string. Some analogies may motivate this assumption. First of all, as is the case with a classical violin string, the Q 
of a resonating string (related to the decay time) can be very large compared to the time it takes to excite ( "pluck" ) 
the string. Secondly, as is the case with cosmic strings, energy loss through particle emission may be inefficient until 
cusps or interconnections are formed. Assumption (5) states that the particle emission rate from the string grows 
nonlinearly with increasing L. We believe this is a reasonable expectation upon nature. E.g., the emission may occur 
as a result of string pinching to form cusps, or crossing to form reconnections, etc. The nonlinear nature of pinching 
and crossing is the reason these phenomena are studied numerically on a lattice, rather than analytically. Also, the 
QCD string is non-abelian, so that interactions among string sections are not linear. In fact, one may wonder whether 
our assumption (3), namely, linear string growth per interaction, is sensible. What really drives the model is that 
the baryonic-monopole cross-section grows rapidly upon interaction, and that the particle emission rate grows even 
more rapidly; linearity of the string growth, adopted in the model here for illustration, is not an essential aspect of 
the dynamics. 

According to eq. (5.1), the geometrical cross-section after the first interaction grows to a very large value: 

Inspired by such a large cross-section, one may ask how many interactions occur within the first-struck nucleus. 
Naively, the answer appears to be extremely large. The number of interactions is given by counting mfp's as the 
monopole travels through the nucleus: 

A^int- / dz[X-'=nNa], (5.13) 

Jo 

where here the appropriate nucleon density is that in a nucleus, nN ^ fm~^. Note that Mnt is a Lorentz invariant, as z 
and A each scale the same way under a boost. For the cross-section above, one gets Ni^t ^ -Estr/A. However, only that 
part of the string within the geometric cross section of the nucleus interacts. Furthermore, there is a causal limit on 
the rate at which the chromomagnetic string expands, and therefore on the rate at which the monopole cross-section 
grows. For a relativistic monopole, one expects the struck string to expand at nearly the speed of light. This leads 
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to a much more sober estimate for the number of interactions per first-struck nucleus, A/jnt ^ A. The number of 
interactions per subsequent struck nucleus is similar, as the cross-section saturates at the nuclear size A^/^fm^. On 
the other hand, the number of struck nuclei quickly becomes enormous. With the monopole traveling at near light 
speed, and the geometric cross-section growing also at near c, essentially all of the nuclei in the forward light-cone 
starting at the site of first interaction are struck, until the baryonic-monopole is quickly brought to rest. 

In passing, we note an alternative possibility for a growing strong cross-section might be to enhance (q- monopole)- 
nueleon hard scattering relative to string micleon scattering, so that q monopole recoil in the first interaction provides 
a large stretch to the string. The recoil energy can be significant, up to providing a stretch SL ^ 7A^^, when 

the condition ikf^ ^ 2Emitii^ ^ (PeV)^ holds. A schematic of q monopole recoil is shown in figure (5d). The simplest 
way to increase the ((7-monopole)-hadron cross-section is to increase the color-magnetic charge of one or more q— 
monopole constituents of the baryonic monopole. Increasing one constituent charge by an order of magnitude increases 
the cross-section by two orders of magnitude. Generally, the values of the constituent charges are model-dependent, 
with larger values typically coming from less attractive models. An added complication with a large [q monopole)- 
nucleon cross-section is that the energy losses of baryonic-monopoles propagating through cosmic media and through 
the magnetic fields of the earth, sun, galaxy, etc. may require a re-analysis. 

The baryonic-monopole's mfp (cropjv)"^ for the first hadronic interaction is an observable and therefore of interest. 
We may estimate the unstretched baryonic-monopole "size" ctq by equating the energy stored in the string with that 
due to the repulsive magnetic force between the g-monopoles. This gives 

Lo ■ 

Thus, the unstretched monopole's string length is 



Mi^o ^ |- . (5.14) 



Lo ~ =i , where g = y -j- - 6. (5.15) 



The lateral size of the string is A ^ so that the geometric cross section is 



/200MeV\ 



2 



£70 ~ 6A-^ ~ 60 mb j . (5.16) 

This is somewhat larger than our order of magnitude estimate A^^, but comparable to the cross section of a high 
energy proton. Accordingly, we expect Xmax for the baryonic-monopole to be similar to that of a proton, or even a 
nucleus. 

We remark on a possible signature discriminating the monopole from a proton. The large number riemit ~ M/ Esti 
of soft individual nuclear interactions underlying the monopole-induced shower event, in contrast to one of a few hard 
interactions from a primary proton, may lead to a larger hadron to electromagnetic ratio in the shower; this in turn 
leads to a higher muon content. 

Finally, we remind the reader that the results of this section are predicated on a very speculative model wherein the 
strong cross-section of the baryonic-monopole grows very rapidly after the first interaction, and the particle emission 
rate from excited strings grows even more rapidly. While these dynamics are plausible. Nature may reject them. Then 
baryonic-monopoles will look much like Z-monopoles in their interactions with matter. 
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VI. SUMMARY AND CONCLUSIONS 



The challenge to seek a cosmic magnetic monopole flux below the Parker limit has now been met technologically 
in various underground/ice experiments. Here we provide theoretical motivation for further efforts to increase the 
search sensitivity. 

The Kibble mechanism predicts possibly observable fluxes for monopolcs with masses in the intermediate range, 
-|^q8±3 QgV. Such monopoles will be naturally accelerated by cosmic magnetic fields to relativistic energies. The Uni- 
verse is transparent to relativistic monopoles meaning that a cosmic monopole flux should arrive at earth unattenuated 
and enter terrestrial detectors. 

The possible signatures by which a relativistic monopole flux may be identified are varied. The Dirac condition for 
the monopole electromagnetic charge is au = 1/4:01. This large coupling makes promising a search for electromagnetic 
signatures including direct Cherenkov emission, and coherent radio-Cherenkov emission from the charged e~-e+ 
shower. The radio signal may be observable for monopoles in polar ice. Of interest is the possibility of "m-raying" 
the interior of the earth with the monopole flux if their mass is within the range 3 — > 10 PeV. 

The hadronic properties of monopoles are not well known. Monopoles with color (g-monopoles) and without color {l- 
monopoles) are expected to be produced in phase transitions which break the gauge symmetry to the Standard Model. 
Just as the monopole magnetic charge is dual to electric charge, the monopole chromomagnetic charge is dual to the 
familiar chromoelectric charge. The Z-monopoles carry QCD fields in their unbroken core, and so have a soft strong- 
interaction of size a ^qcd- '^^^ monopoles may have a more interesting strong-interaction. Presumably they are 
confined by chromomagnetic strings into color-singlet baryonic-monopoles, with initial cross-section ctq ~ ^qcd- 
the process of interacting, the strings may not break until stretched to length at least L ^ M//z, where /i ~ ^qcd 
the string tension. Thus, this stretching (if it can take place) will lead to a much larger cross-section if the energy-loss 
mechanisms of the string are slow compared to the interaction time. We have speculated that this enhanced cross- 
section may be many times larger than the usual QCD cross-section, and if so it could very quickly transfers 0(1) 
of the monopole energy to the shower. In this case baryonic-monopoles would be viable candidates for the observed 
cosmic rays above Egzk ~ 5 x 10^^ eV. 

VII. APPENDIX: SHOWER DEVELOPMENT MODEL 

A monopole is highly penetrating and, as such, can initiate many subsliowers before stopping. However, for a 
subshower to develop, the energy injected into the absorber in any single interaction must be greater than E^- 
Thus, the inelasticity restriction 77 ~ ~ 10~^^ ( io^°ev ) ^ holds for monopole-matter interactions which initiate 
subshowers, with lower inelasticity events contributing only to ionization. 
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Incident Monopole 

\ 



Vacuum 




FIG. 6. A schematic representation of the variables entering into a calculation of the shower development for a 
monopole-induced shower. The monopole path is along the dashed line. The depth where the shower size is evaluated is 
at X. The smallest depth where an initiated subshower contributes to the total shower is a;min. 



Suppose the monopole initiates the j subshower at a depth Xj . The subsequent development of the j -subshovifer 
has particle number 



(7.1) 



where < (x — Xj) < i^e In (^ i^e ) ' total energy injected into the shower at the point Xj is Ej, and A^o is the 
initial number of particles injected into the shower at the point xj (For electron pair production A^o ~ 2, whereas for 
bremsstrahlung and monopole-electron elastic scattering Nq = 1.). 

The total shower development N{x) is then obtained by summing over all the subshowers 



N{x)=J2n,{x). 



(7.2) 



The approximate distance between the initiation of subshowers is given by the monopole mean free path, A = 
where a is the monopole cross-section and n is the number density of scattering centers. The j^^ interaction then will 
roughly occur at the depth Xj — jX. If the inelasticity per interaction is approximately given by a constant value of 
r] then the energy injected at the j^^ scattering is 

V 



1 — T] 



{l-r,yEo. 



(7.3) 



The monopole will continue to initiate subshowers until [77] it is degraded in energy to the point where Ej ~ NqEc- 
That fixes the maximum number of subshowers to be 
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In 



/ (l-r,)NoEA 
\ vEo J 



Note that this relation is just a restatement of the condition < (a; — Xj) < In (^ i^e ) ■ ™ (^•^) 

now be performed from j = 1 to j = jmax- 

The sum in eq. (7.2) can be transformed into an integral over column depth by making a continuum approximation. 
The substitutions X^,- ^ / ^ and xj x' give 



r dx' 



N{x) = No I ^e"T7^ e^^T^ - 1 . (7.5) 



The physical interpretation of the limits of integration is as follows: Xmin is the smallest depth at which an initiated 
subshower is still cascading at the depth x. Subshowers initiated at depths < Xmin will have already ranged out by 
the point x [76] , and so they cannot contribute to the particle number at x and are excluded from the integral. 
From the geometry of fig. (6), using eqs. (4.3) and (7.3) and taking j ^ f > we can deduce 

x^^ix) =x-Xl^^ = x-^,ln j . (7.6) 

Solving for Xmin gives 



^min — c , / \ 

1 + ^ ln(l - r]) 

Using <C 1, so ln(l — J?) — — r? and ~ 77, this result simplifies to 



(7.7) 



X — X 

^min = ^ • (' -oj 

^ A 

Here and in the remainder of this section we define X = Xjjj^x = In (^ nI^^e ) ' which is the depth of the first 
subshower maximum (coming from a monopole with its full kinetic energy Eq). Expression (7.8) is formally singular 
at ?7 = A/^e ^ Ij but for monopoles this singular point is not relevant. 

For a monopole passing through matter, ^ <C 1, so we can write eq. (7.8) as 

X - X^in =X- ^x . (7.9) 

Substituting the above into eq. (7.5) gives the quasi steady-state shower size 

Ar(a;)^|exp(Q (7.10) 

Notice that the a;-dependence has vanished, as is appropriate for a steady-state phenomenon. Using the definition 
of X above, and the useful relation between the continuum and discrete expressions for energy loss 



r]E_ 



dE 



dx 



we rewrite (7.10) as 

dE ■ 



Nil) = ^ 



dx 



(7.11) 



(7.12) 

rtcTit<ri 
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where 7?„it = ^2QMeV_ 




FIG. 7. Electron-positron pair production of an M = 100 TeV relativistic monopole in air where the inelasticity is restricted 

to 77 > TJcrit = ■ 

We next turn to some details of the shower energy. In fig. (7) wc plot the dominant terms for energy loss to 
electron-pair production. X2 and X3- For 7 ~ 10*, pair production is dominated by X2; which describes slow pairs in 
the no screening limit. For shower development we must restrict the inelasticity to r] > Tjcrit = ^oSMeV, xhe plot X2 
shows the effect of this restriction. Notice that X2 — X2 for 7 ~ 10^, so we are justified in using eq. (3.8) without the 
■q restriction to compute the shower size when 7 ^ 10^. In the evaluation of eq. (7.12) we use a parameterization of 

in difi'erent media (good to ~ 2.5%) given by [78] 

716.4 A 2 /„,ox 

cm-2 (7.13) 



and a fit to the data [79] for Ec : 



Z(Z+1) In(^) 



Ec I 7l\'^2d solids and liquids 



MeV 1 for gases 



(7.14) 



Z+1.24 

710 
Z+0.92 

The contribution of the photonuclear process to the electromagnetic shower is indirect. The photonuclear interaction 
injects hadrons into the monopole shower. A subshower initiated by a high energy hadron will produce tt^'s as 
secondaries, which each decay to 2 7's. If these 7's have E > Ec, they may initiate an electromagnetic shower. So, only 
the largest inelasticity fraction of the energy lost via the photonuclear interaction contributes to the electromagnetic 
shower in the end. Given this fact, it is reasonable to assume that pair production alone provides a lower bound to 
the electromagnetic shower size and that the pair production plus photonuclear interaction provides an upper bound. 
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Boost Factor, 7 



FIG. 8. The monopole-induced shower in ice for a monopole of mass 100 TeV. The shower size is the total number of electron, 
positrons, and photons. The dashed line is the total shower from pair production and the photormclcar interaction. The solid 
lines show the electric charge excess (roughly 20% of the shower size) for pair production alone (oc 7) and pair production plus 
photonuclear (oc 7^^'^*). 



3. Lateral Shower Profile 



For monopole-induced showers the lateral profile is greatly simplified in comparison to hadronic primaries. As 
shown previously the monopole-induced shower is approximately constant, being continuously regenerated by the 
monopole as the lower energy particles range out. For our purposes in the following section it is sufficient to assume 
that the lateral profile is uniform out to a lateral cutoff given by the Moliere radius 

p 74 g , ge ,, 100MeV , 

As defined, the Moliere radius is independent of the incident monopole energy, being determined only by the spread 
of low energy particles resulting from multiple Coulomb scattering. Within a distance Rm of the monopole path, 
90 % of the shower particles will be found [78]. 
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